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We present theoretical results for the underscreened Kondo lattice model with localized 5* = 1 spins coupled 
to a conduction band through a Kondo coupling, Jk, and interacting among them ferromagnetically. We use a 
fermionic representation for the spin operators and expand the Hamiltonian in terms of bosonic fields. For large 
values of Jk, we obtain a ferromagnetically ordered solution and a Kondo regime with a Kondo temperature, 
Tk, larger than the Curie temperature, Tc- This finding suggests a scenario for a coexistence of Kondo effect 
and ferromagnetic order In some uranium compounds, like UTe or UCuo.QSb2, this kind of coexistence has 
been experimentally observed: they order ferromagnetically with a Curie temperature of order Tc ~ 100-ft' 
and exhibit a Kondo behavior for T > Tc- The proposed underscreened Kondo lattice model accounts well 
for the coexistence between magnetic order and Kondo behavior and yields to a new "ferromagnetic Doniach 
diagram". 

PACS numbers: 71.27.+a, 75.30.Mb, 75.20.Hr, 75.10.-b 



I. INTRODUCTION 

The Kondo lattice (KL) model is one of the fundamental 
microscopic models for studying the properties of strongly 
correlated electron systems, and a large amount of theoreti- 
cal work was carried out on this problem in recent years (for 
a review, see Ref. 1 ). This model is widely used to describe 
the physics of intermetallic heavy fermion compounds based 
either on rare earths elements, or on actinides^. In heavy 
fermion materials there are two different types of electrons: 
conduction electrons from outer atomic orbitals, and strongly 
correlated electrons from inner /-orbitals, the later ones being 
generally localized. The KL model describes the interaction 
between these two electronic subsystems in the limit when 
/-electrons are completely localized and form a lattice of lo- 
calized spins. 

Historically, KL model has been proposed to account for 
properties of cerium compounds, where a competition be- 
tween Kondo effect and magnetic order has been experi- 
mentally observed. Such competition gives rise to a rich 
phase diagram with various quantum phase transitions. KL 
model has been proven to be an appropriate tool for describ- 
ing these quantum transitions at different values of external 
parameters such as band filling, pressure, magnetic field or 
temperatures^ 

In most cerium compounds, Ce ions are in the localized 
4/^ configuration corresponding to spin S = 1/2. This lo- 
calized spin couples antifeiTomagnetically, via an on-site ex- 
change interaction, J^, to the conduction electron spin den- 
sity. At very low temperatures the localized spin S ~ 1/2 is 
completely screened by the conduction electrons, leading to 
the formation of coherent Kondo spin-singlet state. Besides, 
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the local coupling between /—spins and conduction electrons 
may give rise to a magnetic order through the RKKY interac- 
tion. This interaction is usually added to the KL model as an 
additional inter-site interaction between /— spinJ^El 

The competition between the magnetic order and the Kondo 
effect was first considered by Doniach^ ''. He proposed 
a phase diagram with a quantum phase transition between 
a magnetically ordered phase and a non-magnetic Kondo 
phase. Doniach phase diagram was later extended^ to include 
the short range magnetic correlations that survive inside the 
Kondo phase*^ . From an experimental point of view, the com- 
petition between magnetic order and Kondo effect has been 
observed in many cerium and ytterbium compounds, which 
yields a set of very rich phase diagrams with various quantum 
phase transitions under pressure (see, for example, ref. 6 and 
references therein). 

In this paper, we focus on the physical properties of ura- 
nium compounds. This is another class of heavy fermion sys- 
tems, which show very rich behavior, quite different from 
cerium compounds. It is peculiar of uranium compounds 
that they exhibit numerous coexistence phenomena, the most 
prominent of which is the coexistence of magnetic order with 
Kondo effect*^ - " ' - or t he co existence of the magnetic 
order with superconductivitjD^'^. We will be primarily in- 
terested in the coexistence between ferromagnetic order and 
Kondo behavior, as, to our knowledge, this effect has been 
somewhat overlooked from a theoretical point of view. 

Let us biiefly describe the experimental situation. The 
first experimental evidence of the coexistence between Kondo 
behavior and ferromagnetic order in the dense Kondo com- 
pound UTe has been obtained long time ago^. More re- 
cently, this coe xisten ce has been observed in UCuo,gSb2 ^ 
and UCoQ,5Sb2^^^^- All these systems undergo a ferromag- 
netic ordering at the relatively high Curie temperatures of Tc 
= 102K (UTe), Tc = 113K (t/Cuo.gS'foa) and Tc = 64.5K 
{UCoQ,5Sb2)- Above the ordering temperatures, i.e. in the ex- 
pected paramagnetic region, these materials exhibit a Kondo- 
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like logarithmic decrease of the electrical resistivity, indicat- 
ing a Kondo behavior. This logarithmic variation extends 
down to the ferromagnetic Curie temperature, T^, suggest- 
ing that the Kondo behavior survives inside the ferromagnetic 
phase, implying that the ferromagnetic order and the Kondo 
behavior do coexist. This coexistence, together with the large 
Curie temperatures, are clearly novel features that cannot be 
explained by the standard KL model"*^. 

Therefore, as a minimal model to describe the Kondo- 
ferromagnetism coexistence, we propose an Underscreened 
Kondo lattice model which, we argue, is appropriate to de- 
scribe the 5/^ configuration of uranium ions. 

The underscreened Kondo lattice (UKL) model consists of 
a periodic lattice of magnetic atoms with 5=1 interacting 
with a spin density of conduction electrons via an on-site an- 
tiferromagnetic Kondo coupling. In addition, the localized 
spins at neighboring sites interact ferromagnetically with each 
other. In this case the Kondo effect does not lead to a com- 
plete screening of the localized spins, and the ferromagnetic 
exchange between the (underscreened) spins may indeed lead 
to the formation of ferromagnetic order. 

We warn that the choice of the model for the electronic 
structure of uranium compounds is a question not settled yet. 
Magnetism in these compounds undoubtedly comes from 5/ 
electrons - this has been proven by many experimental ob- 
servations, e.g. by form-factor studies in neutron scattering. 
At the same time, 5/ electron states in uranium compounds 
are in a crossover region between localized and itinerant be- 
havior, and the degree of localization depends strongly on a 
subtle balance between the electronic structure, the effect of 
correlations and crystal field effects. It is often difficult to 
decide, on the basis of the experimental data, between a lo- 
cal Kondo behavior corresponding to a 5/" configuration and 
a mixed-valence situation. One example is provided by ura- 
nium monochalcogenides: US lies closest to the itinerant side 
for the 5f-electrons, USe is in th e middle and UTe is the 
closest to the localized sidJSEan. Recent photoemission ex- 
periments on UTe have been interpreted as favoring itinerant 
magnetism^', but the magnetic moments deduced from mag- 
netic susceptibility experiments in this compound are close 
to the free ion values of uranium, which implies that the 5/ 
electrons are relatively well localized in UTe Moreover, 
the dual nature of the 5/ electrons, assuming two localized 
5f electrons and one delocalized one, has been considered by 
Zwicknagl et al.^^^^who have obtained by band calculations 
a mass enhancement factor in good agreement with experi- 
ment in UPt^ and UPd2Al^ and by Schoenes et al. '''who 
have carefully analyzed the variation of the localization of the 
5f-electrons with concentration and pressure in diluted U S 
and UTe. The electronic structure of uranium and plutonium 
monochalcogenides has been also studied by DMFT calcu- 
lations So, the appropriate description of the electronic 
structure for uranium compounds is a challenging problem 
and depends strongly on the considered system. Here we re- 
straint ourselves to the study of the UKL model applied to 
uranium compounds such as UTe, when the uranium ions are 
relatively well localized and can be correctly described within 
a 5/^ configuration, in which the two 5/-electrons are bound 



into spin S = 1. 

Besides its applicability to the physics of ferromagnetic 
uranium compounds, UKL model is also interesting on its 
own. It is one of theoretical models which can capture the 
physics of the lattice of underscreened magnetic moments in 
a metal, yet it attracted much less attention than the under- 
screened Kondo impurity model, for which there there exist 
various theoretical studies^'', and an exact solution has been 
obtained by the Bethe ansatJ^SSI There has been only few 
studies of the UKL. The UKL model in the form which we 
use here was first proposed in Ref|30l Magnetism and super- 
conductivity in one dimensional UKL model has been studied 
in Ref. 28 The pseudogap formation in the UKL model has 
been studied in a large-N limit in Ref 29 

The crucial step in our analysis of the UKL model lies in 
the choice of the fermionic representation of localized spins 
5=1. We model them by two degenerate /-orbitals with one 
/-electron each. At each site these two electrons are bound 
into 5=1 due to the strong on-site Hund's coupling. Our 
fermionic representation projects out singlet states and satis- 
fies spin algebra in the triplet Hilbert subspace^'. In other 
words, all possible spin transitions that leave the system in 
the triplet subspace, |5 = 1, 5^) are equivalently described in 
terms of auxiliary fermionic operators. 

We show that UKL model exhibits two continuous tran- 
sitions (more precisely, sharp crossovers): the first one, at 
T = Tk, to a non-magnetic Kondo state with a pseudogap 
in the /-electron density of states, and the second one, at 
T — Tc, to a ferromagnetic state. We find that at strong 
Kondo coupling ferromagnetism and Kondo effect do coex- 
ist. We evaluate the Kondo screening, magnetic moments 
of localized spins and conduction electrons for various band- 
fillings in a wide range of coupling parameters and obtain the 
phase diagram with the regions of Kondo-ferromagnetic co- 
existence, non-magnetic Kondo behavior and pure ferromag- 
netism. This phase diagram can be considered as a ferromag- 
netic "Doniach" diagram for the UKL model. 

The paper is organized as follows: In Sec. II, we present the 
model, introduce the fermionic representation for spin opera- 
tors, derive the Green functions and perform a self-consistent 
analysis. The results of the calculations at zero- and finite tem- 
peratures are discussed in Sec.III. In Sec. III.E, we present the 
ferromagnetic "Doniach diagram" for our model. The con- 
clusion contains the discussion of the main results and the 
comparison with the experimental data for the uranium com- 
pounds. 



II. UNDERSCREENED KONDO LATTICE MODEL. 

A. Fermionic representation of localized spins. 

We first introduce and discuss in detail the fermionic repre- 
sentation for localized spins 5 = 1, made out of two fermions 
on degenerate /-orbitals. The two fermions couple into 5=1 
due to the strong on-site Hund's interaction. We introduce 
a fermionic representation in the constrained Hilbert space 
which contains only triplet spin states, dropping out all spin 
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singlet states. 

The projection of the spin singlet states is justified due to 
the correlation nature of /-electrons. The strong Hund's cou- 
pling favors the triplet states with energy Et with respect to 
the singlet states with energy Eg, and to the states with two 
electrons on a single orbital, with even higher energy Ed- 
Then, considering just states with S* = 1 the transformation 
between different Sz (|1, Sz)) projections are unambiguously 
described in terms of fermionic operators aJ^ 

hifhfi^fy2i + /iVu/1t/2t + fhfnfhf2^+ 

/li/iT/lT/2i):|l, 0)^11,0) 
t f ft 

'U/2, 



fljnflj2i : 1 1,-1 



ll,-l> 

i(/iViT/2V2T + /uAt/2t/2t) ■■ Iia> - |1,0) 



(1) 



i(/iVu/2T/2i + fnhiflhi) : |i, -1) ^ |i,0) 
\{fl^h]fl^l2i + /^/ix/2|/2t) : |1:0) ^ |1,1). 

where and f^a are creation and annihilation operators for 
/-electrons, carrying spin and orbital indexes, a and a{a = 
1,2), respectively. In terms of spin operators, the transitions 
represented by Eqs.Q can be expressed as: 

l-5Ml,0)-.|l,0) 
^Sz{Sz-l) : |1,-1)^ |1,-1) 
ISz{Sz + l) : |1,1) ^ |1,1) 

: |1,0) |1,-1) (2) 
:|1,-1)^|1,0) 
1,1) ^ |1,0) 
i5,5+:|l,0)->|l,l) 



Sz ■ 



Note that there is no transition |1,1) ^ |1,— 1). The equiv- 
alence between fermionic and spin representations should be 
interpreted as the equivalence of matrix elements between cor- 
responding states. 



B. Model. 



Our model Hamiltonian is the following: 



(3) 



The first term represents the conduction band with dispersion 
energy e^. We assume that the band has a width 2D and the 
density of states for conduction electrons is 1/2D in the in- 
terval [—D,D\ and zero otherwise. The operators cj^g.(ck(T) 



correspond to delocalized Bloch states with spin cr, while /i 
is the bare electron chemical potential. The second term de- 
scribes the energy of localized levels and E^, can be consid- 
ered as a fictitious chemical potential, i.e. a Lagrange multi- 
plier for auxiliary /-fermions. The actual value of Eo is fixed 
by a local constraint n/ = X^jcra ^^la ~ ^ for the number of 
/-electrons per site. The third term is the antiferromagnetic 
on-site Kondo coupling, with Jk > 0, between localized /- 
spins. Si — 1, and conduction electrons <Ji — 1/2 spins. The 
spin operators ai can be written in terms of the fermionic op- 
erators in a standard way: 



1 (n^l — n^l ). The last term in Eq.(|3|) is a ferromagnetic inter- 



site interaction, Jh < 0, betweenTocalized /-magnetic mo- 
mentJ^H. 

Using the fermionic representation for both /- and con- 
duction electrons spins, we can rewrite the Hamiltonian ([3]) 
in terms of / and c electronic operators. As the expression 
of the Hamiltonian in terms of fermionic operators is rather 
lengthy, we show here in detail just the most relevant term, the 
transverse part of the Kondo coupling ^Jxicr^S^ + a^S^), 
which reads, when expressed in the fermionic representation: 

^JK{crtSi + cr:[S+^ = |Jk(c+Cj|(/+^/j1|/+j^/j2T + 
faifi^]fi2]f^2] + /ix/ju/i^x/»2T + /^x/iiT/^x/i2i)+ 

caCi|(/i||/iix/^i/,2i + /i[t/iu/i^x/j2i+ 

ftl]fil]fa]fi2l + /il|/iu/i^|/i2T) 



(4) 



Now we define the relevant bosonic fields to describe the 
Kondo effect: we introduce the operator A^^ = X]q ^ta — 
X]q '^tcrfia which couples electrons and auxiliary fermions at 
the same site. Using this definition, we rewrite the Eq.d?]) as: 
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1-3+ 



ljK{-Xj^X:fn,-\l^\:ln,- 



And in order to describe the magnetic properties of the sys- 
tem, we introduce the operators of magnetization for both 
/- and c- subsystems: Mi = Sf = \{nl^ — nj^) and 
iTLi — af — \{n1^ — "-ij.)' respectively. 

We restrict our consideration to self-consistent analysis 
(equivalent to slave boson mean-field description), and eval- 
uate all physical quantities in terms of these bosonic fields. 
We then introduce four real order parameters A^r = (A,;cr), 
M — {Mi) and m = {rrii), where (...) denotes the ther- 
mal average. The non-zero values of (AI) and (m) describe 
the magnetic phase with non-zero total magnetization, while 
a non-zero A^ describes the Kondo effect and the formation 
of the heavy-fermion state. The limits of this approach are 
discussed at the end of this subsection. 

Within this mean field (MF) approximation, and using the 
fact that the total number of f-electrons in each sublevel is 
TT-Q = 1, the Hamiltonian ^ is expressed in terms of the four 
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FIG. 1: Schematic plot of the band structure: for each direction of 
spin there are two hybridized and El (k) bands (thick solid lines) 
and one non-hybridized /-level Eq^ (dashed line). The left panel 
corresponds to down-spins and the right one to up-spins. The thin 
dotted line corresponds to the Fermi level. 



order parameters {Xa), {M) and (to) as 

-\JKT.^o.M>^^)K. + h.c.)+ 



(6) 



where 

Ekcr = £k + ^0 



JKN{m){M) - \JhNz{MY 



A, = Jk<7{M) , a = ±1/2 

Aa- = JK(T{m) - \jK{K){Xg 



(7) 



JHza(M) 



being z the number of nearest neighbors. In the presence of 
non-zero magnetization the bands for up- and down- spins in 
both /- and c- subsystems are shifted from each other: pa- 
rameters Act and Ag. define the energy shift for itinerant bands 
and localized /-levels, respectively. The shifts for /- and c- 
bands are in opposite directions because of the antiferromag- 
netic Kondo coupling between them. 

For each direction of spin, the diagonalization of the Hamil- 
tonian (j6]) yields one non-hybridized /-state with energy Eqct 
and two hybridized bands with the quasiparticle energy dis- 
persion given by: 



Ska ± \ (Eoa - Ekff)" 



8ai 



(8) 



where = —\JK{Xg). The ± sign refers to the upper and 
lower hybridized band. The hybridization gap in the energy 
spectrum ~ 2\/2as — \/2 J/fAj is spin-dependent and is 
only present as long as is nonzero. The schematic plot of 
the band structure is presented in Fig|T] 

The energy spectra E"^ (k) depend on a set of external pa- 
rameters such as the Kondo coupling, Jk, and the exchange 



interaction, Jh, the conduction band filling ric, the temper- 
ature T, and a set of internal parameters, /i and £'o, which 
should be calculated self-consistently. 

The present approach is equivalent to other MP theories 
developed for the 5=1/2 KL model, such as saddle-point 
approximation in the path-integral calculation, performed by 
Coleman and Andrei^, or large-A, slave boson saddle-point 
approximation by Burdin et al-'"'. AH these approximate meth- 
ods require some caution. Particularly, MF approximation 
neglects magnetic fluctuations, and the system may possess 
spurious charge fluctuations such that the exact constraint 
/|/t + /J/i = 2 is actually satisfied only on average. 

The validity of MF and slave boson approaches for KL 
model with localized spin 5 = 1/2 has been extensively dis- 
cussed in Refs. I34I35I It has been shown that the MF proce- 
dure becomes exact when the original model is extended to 
N ^ oo different flavors of localized electrons. Still, it was 
found in the same paper that the MF procedure captures the 
correct low-temperature physics even for the actual case of 
N = 2 {S = 1/2)). By the same reasons, the MF approxima- 
tion should work in our case where = 3 (5 = 1). 

Below we will study the model (|6]l by employing equations 
of motion for the Green function (GF) for localized and itin- 
erant spins. All propagators will be written in the fermionic 
representation. 



C. Green functions 

We introduce the following notation for the retarded Green 
functions: = ((/«.; /D), F^^ = ((/a.;/^,)) for 

localized /-electrons ; G^^ = ((cko-; cj^^)) ^'^^ itinerant c- 
electrons, and = ((c^;/L» and G"^^ = ((/a<T;4» 
for the "mixed" cf- states. As it is easier to solve the self- 
consistent system of equations in the momentum space, we 
evaluated all Green functions in k-representation. The equa- 
tions of motion for these four Green functions are given by 



iuj-E,^)F-^{k) = l + agG-,M 
{cj-E,,)F;^Jk)^agG-^{k) 

{cj - £ka)G-Jk) = a^(F-Jk) + F^Jk)) 

{iu - £k.)G^,(k) = 1 + a^(G^,(k) + G^,(k)) 



(9) 



From equations of motion (Eqs.(j9])), by straightforward but 
nevertheless tedious algebra (which we report in the Ap- 
pendix) we obtain the following expressions for the Green 
functions: 
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GccO^) = ~w7(k) 



2W„(k) 
1 



,) 2W„(k) 



.-B-(k) 

£k^-g^(k) 



w-£;^(k) 



-£;^(k) 



-_E-(k) . 
£„„-g^(k) 



(10) 



-B-(k) 

_ £o„-g^(k) 
w-_E^(k) 
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where E^. (k) are given by ^ and 

W^{k) = E^ik) - E^(k) - J (Eoa ~ Eka)^ + 8al .(11) 



The Green functions given by Eqs.( 10 1 will be evaluated self- 
consistently in the next subsection. 



D. Self-consistent equations 

We next construct a close self-consistent scheme to evaluate 
the bosonic fields A|, A|, M and m together with the chemical 
potential f.i and the Lagrange multiplier Eq- This can be done 
by imposing (i) constrains on the total number of /-electrons, 
iif — 2, and c-electrons, ric, correspondingly, (ii) the relation 
between the Fermi surface volume and the number of parti- 
cles, (iii) the relation between the total magnetization and the 
number of electrons, and (iv) self-consistent equations for A| 



and A 



The numbers nj and Uc are expressed via the integrals of 
the imaginary parts of the corresponding Green's functions: 



JE 



(12) 



where the summation is over all k points of the Brillouin zone 
and the orbital indexes a = 1,2. As the /c-dependence in 
Eqs.([8]) and ( [T2| comes only through the bare conduction- 
band energies ei^„, we can substitute the summation in the 
fc-space by the integration over e^. The integration is over 
the interval [— D + Ao-, D + A^] where the density of states 
po — is a constant. 

We first consider the case T = when the calculations 
can be done analytically. Let's first evaluate the number of 
/-electrons. For each polarization of spin, /-electrons can 
occupy one localized level Eq^ and a fraction of the two hy- 
bridized bands i?J(k). The expression for spin-up nj^ can be 
then written as 



2D 



D J-_D+A| Wr(er) 



2D 

7 

4 



iD 



^ 4D J/jf 



(13) 



All parameters in ( [T3| ) were defined in (|7]i and ( [TT] l. We have 
also used the following relations: 



£| - Elie^) = i(£| - Eo^ + W^t(^t)) 
£T-4(£t) = K^T-^ot-W^t(£t)) 

Taking into account that f^f" = Jx^ + c 



(14) 



we obtain: 



T 



AD 



1 
AD 



{D + ^^-Eo^Y 



y(^t+A|-£;oT)2 + 



(15) 



The quantity present in Eqs.( T3p5 i as the limit of Integra 



tion over e|, is related to ^ through the relations = /i. 

By the same procedure we obtain the expressions for n|, 

il and nl: 





Mi 

AD 


1 

AD 






-A^-i?oi)2 + 8a2 


(/"i - 






+ 8a^ 


1 


= ^ 4 


Mt 
" AD 


1 

AD 


y 


{D + - Eo^)^ + Sal 




F A| - 




+ 8al 


1 



"■c A 



AD 



1 

AD 



^(Aii + A^-i?oi)2 + 8a 



(16) 



where /i| is obtained from = ji. 

Now we can construct the system of self-consistent equa- 
tions. The first one is the constraint on the total number of / 
electrons: 



n 



7 - 

As we already discussed above, we replaced the local con- 
straint rifi = 2 at each site i by a softer one, for the average 

Uf. 

The second self-consistent equation is obtained by setting 
the average number of conduction electrons to be equal to the 
filling value n^- 

Tic = n[ + n\. . (18) 

The third and forth self-consistent equations are obtained 
from the Luttinger theorem, and from the condition that the 
total magnetization Mtot is the sum of average magnetization 
of /-and c-electrons. This gives 



(17) 



and 



M + m=l 



MT +m 
2D 



- Mi 
4D 



(19) 



(20) 



We used here M = \{n^f — n^j) and m = \{n]. — n^), re- 
spectively. We remark that the Fermi surface encloses both 
conduction electrons and partially localized /-levels. 

The last two equations are self-consistent relations for the 
bosonic fields A| , A^ . These two fields are related to the imag- 
inary part of the mixed c/-Green function: 

A, - /^_^^_^ dujY.^[^l{ImG-,{^) + /mG-2(k))] , (21) 

Using ( [T0| we then obtain the following expressions for A^: 



At 



O "I 
~^2D 



D-\-Ai 



der 



D+A^ W^{e-f) J~D+Aj W^(e^)i 

D+A-f de-f ('22) 

D+Ai Wl{6^) 



— _2LL r 

D J/if+Ai W^{e^) 
oij rMi+A; dei 
D J- 
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Using J 



dx 



ln(a; + \/2? + c) we find 



-D+Af -_Eof +H'f (-D+Af ) 
p-r+AT--EoT+M'T(MT+^T) 

At^+A|-go|+W|(Mi+Ax) 
-r>+Ai-_Eoi+Wi(-D+AJ 



(23) 



Equations ([T7]i,([T8]l,([T9|l,(|20jl and Q constitute the 
full set of self-consistent equations for six variables: 
A|, A|, M, m, and Eq. We solved this equations numer- 
ically by iteration and we explicitly verified that the solution 
minimizes the total internal energy Etot 



1 

4D 



Mt-^ot + (Mi + D)Eqi- 



1 

8r> 



{D + A^-Eo^)W^{D + A^)- 
{fi^ + A| - £'o|)VF|(/i| + A|) 
+ p + At + i?oi)W^i(-i? + A^)+ 



(24) 



0^4+A4-£;oi)M^i(A/i + A4) 



-JKmM - IJhzKP - EpUc - 2Eo 



At finite temperature, T, the set of self-consistent equations 
remains the same, however we are no longer able to obtain an 
analytical expression for the occupation numbers and bosonic 
fields. The relations between Uj, n^, Acr and the Green's func- 
tions are now: 



n 



< = d^npicu) Ek[-i/mG-,(k)] (25) 
A. = /+^ diunpiiu) E^i-Hl^G^M + /mG-2(k))], 



where np{uj) 



e T +1 



is the Fermi distribution function. 



Straightforward calculations lead to 



n 



£o<,-S"(£<,)i 



(26) 



+nF{E'L{e„)) ^^^^^^ , 

A- = S S^'otL dea[nF{El{e,)) - np{E-_{eMw^) 
As for T = 0, we now evaluate numerically the integrals 



in Eqs.( 26 1 and verified that the numerical solution of self- 



consistent equations minimizes the free energy of the system 
F. The latter can be calculated through the partition function 
Z: 

Z = ^[e-/3(^-(k)-M) + e-^(^?('^)-'^)] (27) 




FIG. 2: Plot of the parameter A| as a function of Jk and Uc for 
Jh = —0.01. There is a discontinuous transition as a function of 

1/2 

Jk and a behavior in ric as explained in the text 



and 



F = -TlnZ 



(28) 



In the following section, we present the detailed discussion 
of the results obtained with this self-consistent scheme for the 
UKL model. 



in. RESULTS. 
A. T=0 

We first discuss the properties of the model at T = 0. In 
order to establish the region of coexistence of Kondo effect 
and ferromagnetic ordering (Kondo-ferromagnetism coexis- 
tence) in the phase space set by the parameters of the model, 
we study the behavior of Acr and the magnetization M and 
TO. These order parameters are correlated and it is useful to 
discuss them together In Fig|2] we present the variation of 
the Kondo correlation A| as a function of the Kondo coupling, 
Jk, and the band-filling, ric, for a fixed value of Jh — —0.01. 
In FigjSjwe present the dependence of the total magnetization 
(A/ + to) on the same parameters. In these and all other fig- 
ures all energies and temperatures are measured in units of the 
half-bandwidth D. The values of A^ and A| are close, there- 
fore we present just A| in Fig. |2] in order to keep the figure 
transparent. 

From the analysis of both figures one can verify that for 
values of Jk smaller than 0.6, the Kondo correlation A| re- 
mains equal to zero (Fig|2]) for all values of Uc considered. 
The total magnetization Af + m is large and equal to its max- 
imum value (FigjSji. This phase is a pure magnetic one with 
no Kondo effect. When Jk increases, for a fixed value of ric, 
there exists a critical value of Jk above which both A| and 
the total magnetization Mtot change abruptly and A| become 
different from zero. The region of the parameters Jk and ric 



FIG. 3: Plot of the total magnetization {M + m) as a function of Jk 
and ric for Jh = —0.01. There is a discontinuous transition as a 
function of Jk that corresponds to the transition to A| = in Fig[2] 



where A| and the total magnetization Mf at are finite is the re- 
gion of coexistence between Kondo effect and ferromagnetic 
order. As the magnetization of localized and conduction elec- 
trons, M and m, have opposite signs due to antiferromagnetic 
Kondo coupling between them, the maximum value of the to- 
tal magnetization is always less than 1. The critical value of 
Jk decreases with increasing ric, as it is expected from the 
"exhaustion principle'f^SEl. 

One can see in Fig|2]that decreases smoothly as a func- 
tion of ric, while it undergoes a sharp transition as a function 
of Jk. We calculated the variation of A with and obtained 
an approximate square root behavior A ^ {ric — ri c.crY, with 
r w 1/2, similar to the KL case with spin 1/^51211, Indeed, in 
Fig. [4]we present a plot of A| as a function of Uc for Jk — 0.9 
and Jff = —0.01. The solid line in this figure is the best fit to 
the numerical results and yields a power law behavior with an 
exponent very near 1/2. 

We next study the effect of the exchange interaction be- 
tween localized spins, Jh, on the region of the Kondo- 
ferromagnetism coexistence. In Fig|5]we present the phase 
diagram in the {Jk, Jh} plane for a fixed value of ric- The 
region of coexistence extends up to Jk — 1-2 and down to 
Jk = 0.5 while Jh can vary between —0.001 and —0.06. 
The values of Jh in Fig |5] are rather small; however we recall 
that the real strength of the "local field" applied on an individ- 
ual spin by its neighbors is zJh, where z = 6 is the number 
of nearest neighbors in a simple cubic lattice. 

Summarizing, at T = there exists a region of coexis- 
tence between the Kondo effect and the feiTomagnetic order 
in a wide range of parameters Jk, Jh, and Uc, including the 
case of the half-filled conduction band, ric — 1. Also, a dis- 
continuous transition to pure feiTomagnetic state is found by 
decreasing Jk while a continuous change is obtained varying 
the band filling, ric. 
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n 



FIG. 4: Plot of the parameter A^ as a function of Uc, for Jk ~ 0.9 
and Jh ~ —0.01. The circles correspond to the numerical calcula- 
tion and the line to a fit with a square root law. A| goes to zero for 
He ^ 0.42 
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FIG. 5: The dark region indicates the values of the parameters Jk 
and Jh where there is coexistence between the magnetic ordered 
state and the strongly correlated (Kondo) state. The intensity of the 
gray tone is proportional to the value of A| 



B. Finite temperatures 

We now present the results obtained at finite temperatures. 
First we show in Fig.|6]the temperature variation of the Kondo 
correlations A| and A|, and also the / and c magnetization M 
and TO. The parameters used in the calculation of Fig. |6] are: 
Jk — 0.8, Jh = —0.01 and ji^ — 0.8. The two magnetiza- 
tion curves clearly show the existence of a continuous phase 
transition at the Cuiie, Tq, from a ferromagnetic to a param- 
agnetic state. At low temperatures we observe the coexistence 
of a magnetic order and Kondo behavior The strength of the 
Kondo effect is the highest at Tc- As a matter of fact, when the 
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FIG. 6: Plot of A|, A|, M and m as a function of temperature for 
Jk = 0.8, Jh ~ —0.01, ric — 0.8. At zero and low temperatures 
we observe the coexistence of magnetic order and heavy fermion be- 
havior, the Kondo effect is maximal when magnetization vanishes. 



magnetization decreases, grows and passes through a max- 
imum at the Tc; its value at that temperature is about 20% big- 
ger than at T = 0. Due to the breakdown of the spin symme- 
try, and A| are slightly different in the magnetic region but 
they coincide at Tc, when the spin symmetry is restored. For 
T > Tc the system exhibits only Kondo behavior (A^ ^ 0, 
M = and m = 0). Finally, we define the temperature at 
which Act vanishes and / and c electrons become decoupled, 
as the Kondo temperature, Tk- The fact that A vanishes at a 
particular temperature, instead of slowly decaying to zero, is 
a well known artifact of the mean-field approximation. Actu- 
ally, Tfj- is a crossover temperature, associated with the onset 
of local Kondo screening. 



C. Density of states 

We also analyze the /- and c- electron densities of states 
(DOS) at various temperatures. The /- and c- electron DOS 
are calculated numerically from the imaginary part of the /— / 
and c — c Green functions. 



A 7^ 

koL 



(29) 
(30) 



Fig|7] shows the / and c DOS at T = 0. The parameters 
are the same as in Figure [6] For these parameters, both the 
magnetic order and the Kondo effect are present in the ground 
state. We see from the figure that the bands for the two spin 
polarizations are shifted, as expected for a magnetically or- 
dered state. The hybridization gap due to Kondo effect is 
present for both directions of the spin, but they do not coincide 
neither in the borders nor in the width. The non-hybridized /- 
level lies inside this gap. This localized /-level is occupied 



FIG. 7: Schematic plot of the density of states for the magnetic 
case at T = 0. The parameters are Jk = 0.8, Jh = —0.01 and 
Tie — 0.8. The continuous line corresponds to the /-density of states, 
the dashed line to the conduction density of states and the dotted ver- 
tical line indicates the position of the Fermi energy. The peak of the 
localized f-level lies inside the gap and here the up localized f-level 
is occupied and the down localized f-level is empty. The inset shows 
the variation of the hybridization gap for the two spin directions as a 
function of the temperature. 



for spin-up states and it is empty for spin-down states. We 
caution that the finite width for the two localized /— levels 
in the figure is due to the fact that we added a very small but 
finite imaginary part to the energies. The inset of Fig|7] shows 
the temperature variations of the hybridization gaps for two 
spin directions. The two gaps are nearly equal because they 
are proportional to X^, and the difference between A| and A| 
is rather small. 

In Fig|8]we plot the density of states at four different tem- 
peratures: a) r = and b)a low-temperature phase {T < Tc) 
both exhibiting coexistence of magnetic order and Kondo be- 
havior, c) pure Kondo phase at Tc < T < Tk and d) un- 
corrected high-temperature phase at T > Tk- In the coex- 
istence region at T < Tc, the Fermi level remains inside the 
hybridization gap for the spin-up band, and inside the con- 
duction band, E'^{k), for the spin-down band. That implies 
a semi-metallic behavior. When T > Tc, the up- and down- 
spin bands coincide, the Fermi level lies inside the gap and 
coincides with the energy of the /-level, Eoa- That implies 
an insulating behavior Finally, when T > Tk, Xa vanishes 
and there is no more coupling between / and c electrons: the 
hybridization gap closes and the system becomes metallic. 



D. Effective mass. 

^From the quasiparticle spectrum ([8]l one can estimate the 
mass enhancemenP^ 



m 



2a 



(31) 
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FIG. 8: Plot of the up- and down- f-density of states for T = and 
three finite temperatures 
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FIG. 10: Plot of the electronic specific heat for the parameters a) 
Jk = 0.8, Jh = -0.01, Uc = 0.8 (circles) and b) Jk = 0.7, 
Jh = -0.001 and Uc = 0.8 (triangles). 
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FIG. 9: Plot of the effective mass versus temperature with the same 
parameters as in Fig|6] The triangles correspond to up- spin and the 
circles to down- spin. 



We plot the temperature variation of the mass enhancement 
in Fig|9] Two peaks are clearly seen: one corresponds to the 
Curie temperature, and the second (at a higher temperature) 
to the onset of Kondo effect. In the region of coexistence, the 
effective mass increases as a function of temperature track- 
ing the position of the Fermi level inside the E'^{k) band. In 
the pure Kondo phase at Tc < T < Tk, the mass enhance- 
ment becomes quite large. Such a large enhancement is the 
consequence of the fact that the denominator in Eq. ( [3T] i goes 
to zero in the pure Kondo regime, and the effective mass for- 
mally diverges. A finite width of the /—level eliminates the 
divergence, but still yields a very large mass. Although this 
estimate of the mass probably overestimates the strength of 
the enhancement, it still provides a qualitatively good expla- 
nation of the heavy fermion behavior 



E. Specific Heat 

The thermodynamic properties of the UKL model can be 
easily calculated from the expression for the free energy 
Eq. 



(28 I 



Here, we present only results for the specific heat 
which could be compared with available experimental data. 
We use the usual definition of the entropy of the system 



a pi r) 
and the specific heat can be calculated as 



Cv{T) = T 



dT 
2^2 



(32) 



(33) 



1 (£;^(k)-/i)^e^(^-('-)-^) 
t^ , „/3(£;^(k)-M)\2 



1 



(1 + ^l3(.El(\L)-t,)Y 



E 



{Eoa -life 



(1 + 6/3(^^0- -A')) 2 



In Fig. 10 we present numerical results for the specific heat 
as a function of temperature. As the specific heat is very sen- 
sitive to the choice of the parameters, we plot it for two dif- 
ferent sets: a) Jk — 0.8, Jh = —0.01, ric = 0.8 (circles), 
and b) for Jk = 0.7, Jh = -0.001, n, = 0.8 (triangles). 
For the later choice of parameters, the specific heat tempera- 
ture dependence is more pronounced. First, at low tempera- 
tures, the specific heat increases with increasing T, following 
a power law with an exponent close to 2. Second, it goes 
through a peak, which can be associated with the temperature 
when both A| and A| become equal and to the Curie temper- 
ature. At higher temperature, in the non-magnetic region, the 
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FIG. 11: The ferromagnetic Doniach diagram : Plot of the Curie 
temperature Tc (the triangles are analytical values, the circles nu- 
merical ones) and the Kondo temperature Tk (squares) versus Jk 
with Jh = -0.01 and = 0.8. 



gram (localized spins S=l/2), the reduction of the Neel tem- 
perature at large Jk is the consequence of the competition 
between magnetic order and Kondo effect. The ordering tem- 
perature goes to zero at high enough values of Jk ■ Here, there 
is no such competition, and the Kondo interaction actually fa- 
vors a ferromagnetic ordering through the RKKY interaction 
(even if the localized magnetic moments are partially screened 
by the conduction electrons). That explains the augmentation 
of the Curie temperature. Then, if in the standard KL model 
the short-range antiferromagnetic correlations between neigh- 
boring /-moments reinforce the Kondo effect^ ^, in the UKL 
model it is the Kondo effect that reinforces the ferromagnetic 
interaction. 

The diagram presented in Fig.[TT|can be called the "ferro- 
magnetic Doniach diagram" for the UKL model and it is qual- 
itatively very different from the well-known Doniach diagram 
derived for the Kondo lattice with S — \/2, where the order- 
ing temperature tends to zero at a finite Jk and the magnetic 
order and the Kondo effect compete rather than co-exist. 



specific heat is very high because the Fermi level coincides 
with the localized /-level, as it was discussed for the effective 
mass. Those values are so high that they mask a second peak, 
expected for the Kondo temperature. Therefore, when the gap 
in the density of states closes at the Kondo temperature, only 
a change in the slope is observed. In the case a) both peaks 
are "rounded up" because of the higher values of the order 
parameters compared to the case b). 



F. Ferromagnetic Doniach diagram 



We present here the ferromagnetic "Doniach diagram" for 



the UKL model. In Fig. 11 we plot the Curie temperature, 
Tc, and the Kondo temperature, Tk, versus Jk for fixed val- 
ues of Jh and ric- It is possible to see that the Kondo temper- 
ature Tk becomes finite only at the critical value ~ 0.65 
for Jh = —0.01 and Uc = 0.8, then rapidly increases for 
larger values of Jk- On the other hand, the Curie tempera- 
ture, Tc, is finite for all studied values of Jk- The two curves 
Tk{Jk) and Tc{Jk) cross slightly above and for larger 
values of Jk the Kondo temperature, Tk, is always larger 
than Tc- Indeed, the ferromagnetic order persists for all val- 
ues of the ratio Jk / Jh, while the Kondo-ferromagnetism co- 
existence exists only for sufficiently large values of this ratio. 

In the purely magnetic region the Curie temperature can 
be easily evaluated, and it is equal to Tc — Also, 
the Kondo temperature exhibits an almost linear behavior as a 
function of Jk- This is in contrast with the Kondo impurity 
case, where an exponential dependence on Jk is observecPH, 
and also with the 5=1/2 KL model where both exponential 
and non-exponential behaviors can be obtained^. 

It is worth to note that the present results indicate that the 
Curie temperature increases as a function of Jk, and this is 
also in the opposite direction of the standard Doniach's dia- 
gram. This can be understood because in the Doniach's dia- 



IV. CONCLUSIONS. 

In conclusion, in this work we introduced and studied the 
underscreened Kondo lattice model with localized spins S = 
1. In the framework of this model, we analyzed the coexis- 
tence of ferromagnetism and Kondo behavior. To this end, we 
derived and solved a set of self-consistent equations for the 
relevant bosonic fields. We obtained a region in the space of 
parameters, where the order parameters, A^, M, and m are si- 
multaneously different from zero. This coexistence persists up 
to finite temperatures, then the magnetic order disappears first, 
at the temperature Tc- At this temperature Acr is still different 
from zero. In fact it is at its maximum value, indicating that 
the Kondo effect is the strongest. This is, in fact, in agreement 
with experiments for uranium compounds in which the Kondo 
behavior is observed above the Curie temperature. The phase 
transition at Tc is, therefore, a transition from the ferromag- 
netic Kondo phase at low T to the non-magnetic Kondo phase 
at higher T. Another transition can be obtained at T = by 
varying the Kondo coupling parameter Jk- This transition is 
discontinuous and leads to a non-Kondo magnetically ordered 
state (see Fig.|2]i. 

We emphasize that the coexistence of ferromagnetism and 
Kondo effect is due to the presence of two /—levels: one 
non-hybridized and one hybridized levels (see Fig. |7]). The 
presence of the non-hybridized /—level is at the origin of the 
partial screening of the localized spins. Indeed, this is the 
key difference between our UKL model and 5 = 1/2 Kondo 
lattice model. In the latter, there is only one /—level and it 
is always hybridized with the conduction band in the Kondo 
regime. The result is a competition between magnetism and 
Kondo effect which is reflected in Doniach diagram. 

Our UKL model can explain the behavior of some ura- 
nium compounds such as the previously described UTe, 
UCuQS)Sh2 and UCoQ,5Sb2 compounds, which order ferro- 
magnetically at a large Curie temperature, Tc, and present a 
Kondo behavior Such a Kondo-ferromagnetism coexistence 
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has been recently observed in other uranium compounds. 
We would like to men tion UN 1812 with a Curie tempera- 
ture Tc = 95iil2MMU, UCoo,QNio,iSi2 with Tc = Q2I^ 
and U Ru2-xRcxSi2 compounds where Tc increases rapidly 
with concentration x^^ In these compounds Tc increases 
rapidly with x and there are clear evidences of the coex- 
istence between the ferromagnetic order and a Non-Fermi- 
Liquid behaviour"*^. A more extended analysis of the UKL 
model can, possibly, describe also this phenomenum. Our 
UKL model can finally be applied to the case of the re- 
cently observed neptunium compound NpNiSi2, which be- 
comes ferromagnetic at Tc = 51.5K and presents Kondo be- 
haviour"*^. In the neptunium based compounds, 5/ electrons 
are relatively well localized and the magnetic moment of Np 
can be described by a localized spin larger than S — 1/2, 
again corresponding to underscreened case. 

In conclusion, the UKL is a very good tool to describe some 
ferromagnetic-Kondo compounds and can be considered as an 
improvement with respect to the regular Kondo lattice model 
to account for some actinide 5/ compounds. 
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VI. APPENDIX 

In this Appendix we present the derivation of the Green 
functions (Eqs. ( 10 1) using the equations of motion (Eqs.(j9]l). 
We first derive the Green functions for the f-electrons: 



(J^.;o(k) + F^„(k)) 



(34) 



where we call Ag- 



Then 



,Eq. Q 



becomes 



(35) 



Solving this set of equations, we obtain 



(36) 



The poles of F^^{k) determine the spectrum E^{k), given 
by eqs. ([sj) in the main text. In terms of i?J(k), the Green 
functions of /—electrons are: 



1 



(i^-Eo„)(Lo-Eo„-2Ag) 

1 rek^--E^(k) £k^-B^(k) 



(37) 



ij-EZW 



and 



1 1 



£k.-S^(k) e^^-EICk), (38) 



2{u:-Eq^) 2W„(k) L w-£;-(k) 



-S^(k) 



where W^(k) ee £:^(k) - E'L{k). 

We next calculate G^„(k). Using the corresponding equa- 
tion of motion 

- ek.)G-Jk) = a^(F-Jk) + (k)) , (39) 



we obtain: 



(t^-ek<T)(i^-Bo<T-2As) 
1 1 1 



(40) 



W<,(k) Lw-_B"(k) w-_E^(k)J 



Finally, the Green functions for itinerant electrons G%^{k) 
are calculated in the same way, starting with: 



{lo - ek.)G^,(k) = a,{GU^) + G^,(k)) + 1 
{uo - Ea.)Gl,{k) = a-,GZ,{k) 



(41) 



And solving eqs. (41 1, we obtain 



rBo^-£^(k) _Eo„-S"(k), 



ui-EZik) 



(42) 
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Nevertheless, we should notice that a kind of coexistence of fer- 
romagnetism and Kondo-like behavior has been observed too in a 
few Ce-based dense Kondo systems, such as CePt^^Si compouncP^ 
or CeAg one". But the Curie temperatures of these compounds 
are relatively small, typically of order 5A", and this result can be 
considered as a clear sign of a strong competition rather than a real 
coexistence between the Kondo effect and ferromagnetic order 



